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Abstract. We seek to obtain a new class of exact solutions of regular black holes in /(T) Gravity 
with non-linear electrodynamics material content, with spherical symmetry in 40. The equations of 
motion provide the regaining of various solutions of General Relativity, as a particular case where the 
function f{T) = T. We developed a powerful method for finding exact solutions, where we get the 
first new class of regular black holes solutions in the /(T) Theory, where all the geometries scalars 
disappear at the origin of the radial coordinate and are finite everywhere, as well as a new class of 
singular black holes. 
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1 Introduction 

General relativity (GR) is formulated in such a way to describe the gravitational effects, that are said, 
the principle of equivalence, equivalent to the inertial by the Riemannian Differential Geometry. The 
curvature is a primary concept of this type of geometry, which should be essentially responsible for 
gravitational effects. The GR has been a theory very effective in describing many local phenomena such 
as black holes and others which are astro-physical, as well as global phenomena, and the evolution 
of our universe. But in particular, the description of the evolution of our universe, by the GR, 
is normally done by the introduction of an exotic fluid, commonly called dark energy, whose the 
pressure has to be negative. This fluid would be the main content that leads to the accelerated 
evolution of our current universe. After several models that attempt to describe the Lagrangian 
associated with the dark energy density, such as the vast majority through the scalar field, there 
are numerous models that, instead of using the GR for this description, modifies it for this purpose. 
The first attempts were the use of a Lagrangian density that contained non-linear terms in the 
curvature scalar R, and so, through the generalization of this idea, there is the best known and 
addressed modification, which is f{R) Gravity [1]. In this theory the action of the system is given by 
= f d‘^x-^—glf(R) + 2K^£jnatter], where f(R) is an analytic function of the curvature scalar. 
Although there were several other modifications, in which the best known are: the Theory f{R,T) 
[2], with T being the trace of the energy-momentum tensor; the f{G) Gravity [3], with Q being Gauss- 
Bonnet term; the f{R, G) [4]. All these theories are commonly generalizations of GR in the sense that 
at some limits of the parameters, they fall in GR. 

Another possible approach to gravitational effects, is one that uses a description of these effects 
directly by inertia. A geometry also alternative to the Riemannian one, would consider an identically 
zero curvature, but at the same time a non-null torsion of the spacetime. Now using the torsion of the 
space-time as primordial concept to describe the effects of inertia, or equivalently, the gravitational, we 
can formulate a theory completely equivalent to GR. This is commonly called Weitzenbock spacetime, 
where the equations of motion are dynamically equivalent of the ones of GR, with a total divergence 
term more in action. The theory described here is known as Teleparallel Theory (TT) [5]. We can 
construct an additional scalar similar to the curvature scalar for the TT, which is represented by 
T, called the torsion scalar. As the equations of TT lead to the same as that of GR, the problem 
of accelerated evolution of our universe remains the same. Then there is the possibility of also 
generalizing the TT, through an action that contains non-linear terms of torsion scalar. The f(T) 
Gravity [6] is the generalization of TT, where f{T) is an analytic function of the torsion scalar. 

Several applications in cosmology were conceived within this f(T) theory, where we can list 
the main advances in [6]. In relation to local phenomena such as black holes, there are not yet 
great advances. The first solution is obtained by Wang in [7], with a coupling with the Maxwell term. 
Afterword, charged solutions also have been obtained in (1-|-2)ZI [8]. Still in the same way, Capozziello 
et al. [9] get charged solutions in d-Dimensions. Vacuum Solutions, known in GR, were re-obtained 
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in [10] with spherical symmetry, and [11] with axial symmetry. The solution of Reissner-Nordstrom- 
AdS (and their particular cases) was re-obtained in [12], and also with new charged solutions in f{T) 
Gravity, where f{T) was reconstructed by a parametric method. 

The initial idea of generalising the TT was the use of an analogy to the action of Born-Infeld 

[13] , where there are non-linear terms for the scalar constructed from the Maxwell tensor, F = 

where F^^" are the Maxwell tensor components. But the original idea of Born-Infeld 

[14] was getting a charged solution where the electric field or by duality the magnetic field, were 
regular at the origin of the radial coordinate. In the Maxwell’s Classical Electromagnetism, the 
solution to a specific electric punctual charge is divergent where this charge is situated, which for the 
description with spherical symmetry, is the radial coordinate origin. This was regularised by the Born- 
Infeld solution. Thus arises the possibility of solutions, coupled with the Non-Linear Electrodynamics 
(NED) Lagrangian density, where there is no singularity in spacetime, which has been noticed in the 
work of Hoffmann and Infeld [15] still at 1937, where they found a regular solution. Long time later, 
Peres obtained another solution where the space-time is free of any type of singularity [17]. A few 
years later, Bardeen suggests a magnetic solution where the space-time is completely regular [18]. 
Thus, new regular solutions emerged [19] and are still fully explored nowadays [20]. 

Exact solutions in GR commonly are obtained containing a singularity, or singular region, which 
prevents the classification of space-time as globally hyperbolic, thus appearing point or non-physical 
region for a given measurement singularity. As in the singularities of the space-time geodesics can not 
be extended, through a finite amount of proper time of a certain measuring apparatus, it disappears, 
which can not be physically acceptable. For this reason, since the beginning of RG, Einstein could not 
accept this kind of singular solution. Until now, it seems that this should only be a classic problem, 
being regularized through quantum effects of gravitation. But one main idea in an attempt to obtain 
regular exact solutions, is that of Sakharov [21], where the imposition “pressure = -energy density" 
(further facilitate obtaining a accelerated current phase of our universe) leads to a breaking at least 
one energy condition, thus satisfying almost all the singularity theorems for that space-time is always 
regular. The motivation to use a Lagrangian density of the NED is arising out from the fact of a 
possible analogy between the effects of polarization of the quantum vacuum, thereby associating the 
probability of Schwinger formula with the component (0 — 0) for the energy-momentum tensor, thus 
being facilitated to obtain a new solution that does not impose a priori the functional form of the 
Lagrangian density of matter [22]. 

Recently new charged black hole solutions have been obtained in /(T) theory, coupling matter 
content from NED, demonstrating the possibility of obtaining regular solutions [23]. Another exact 
solution in scalar-torsion gravity arising of coupling with a scalar field is obtained in [24] for a new 
wormhole-like solution dressed with a regular scalar field. In this work we focus our attention on 
obtaining new regular black hole solutions in the context of f{T) coupled to NED. 

The paper is organized as follows. In section 2 the basis of /(T) Gravity coupled with NED are 
presented, and the equation of motion. In section 3 we show how to reobtain the equations of motion 
of the cases of regulars black holes, arising from the NED, in GR. The section 4 is devoted to the 
presentation of the method of solving the equations of motion arising from the subsections 4.1 for the 
class of solutions with the fix functional form of a(r), and 4.2 for the class of solutions to the form of 
a fix b{r). We make our final considerations in section 5. 

2 The equations of motion in /(T) Gravity 

A differentiable manifold has a structure in which the essential elements of a theory of gravitation 
can be derived from its differential geometry. Through the vectors which are defined in the tangent 
space to the manifold, and co-vectors, dual elements to vectors, which are defined in the co-tangent 

space to the manifold, we can then set the line element of the manifold as 

dS'^ = g^ydx^dx'' = riabe°'e^ , ( 2 . 1 ) 

e“ = e^'^dx^ , dx^ = ej^e'^ , (2.2) 
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where g^^i, is the metric of the space-time, rjab the Minkowski metric, e“ the tetrads ( that belongs 
to the co-tangent space) and and their inverses the tetrads matrices satisfying the relations 
and = <^& • In terms of these matrix, the root of the determinant of the metric is 

given by fy—g = det[e'^^] = e. Then, one can define a geometry where the curvature should be set 
equal to zero. Thus, we have to define a connection for which all the components of the Riemann 
tensor are identically null. This can be done through the Weitzenbock connection. 


pa _ „ a Q 

IJ.U ~ 


„a _ _ a O p o 


(2.3) 


Now we can define two other important tensors to describe the geometry, that of the torsion and 
the contorsion 


rpCX _ -pQ; pQ! 

^ ^ ^ fiy 

(rpliU 

ot 2 ^ 

Through these tensors we can construct a scalar analogous to the scalar curvature of the GR, but 
it is still best to first define a new tensor, from the torsion and contorsion, to simplify the expression 
of the torsion scalar. Thus, we define 


= e,°‘d^e\) , (2.4) 

- . (2.5) 


q ^y _ _ 
- 2 



( 2 . 6 ) 


From the components of the torsion scalar and the ones of the torsion , we finally define the 
torsion scalar as 


rp rpOi q fJ,y 

^ ^ fiy^a. 


(2.7) 


We now obtain the equations of motion taking a mater content as a non-linear electrodynamics. For 
this theory, the Lagrangian density is given by 


C = e[f{T) + 2K^£NED{F)] , 


( 2 . 8 ) 


where Cned{F) is the contribution of the non-linear electrodynamics (NED), and the scalar F = 
{1 /for F^n being the Maxwell tensor components, and = SttG/c^, where G is the Newton 
constant and c the speed of light. The field equations for this theory are obtained through the Euler- 
Lagrange equations. Thus, one has first to take the following derivatives with respect to the tetrads 


dC 

de°-,, 


dry 


= f{T)ee/ + e/T(r)4e“„n„5/‘^ + 2«' 


:dCNED 

de°-,, 


dC 


d{dore°- ) 


= -4/t(T) 5„ (ee/^,'^") - d^TfTriT) + 2^dyy 


dCNED 


didcyF^ ) 


(2.9) 


,( 2 . 10 ) 


where we define friT) = df(T)/dT e fTriT) = d?f{T)/dT'^. Through the above expressions and the 
ones of Euler-Lagrange 


dC 

de 


- da 


A* 

and by multiplying all by the factor e“^e“^/4, we get the following equations of motion 

Sp^‘^daTfTT{T) + [e-^e%da (ee/5,^“) + /t(T) + i^^/(T) = ‘ 

where 7^ ^ is the energy-momentum tensor of non-linear electrodynamics source 


d£ 


_d{dae‘^ )_ 


= 0 , 


n- __ _ 

f R — o 


5^p£ned{F) - 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 
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The choice of the set of tetrads occurs such that all equations are consistent. In our particular case 
we obtain symmetric and static spherically solutions, so we must choose consistently tetrads with 
this final order. If we choose a set of tetrads diagonal as [e“^] = r, r sin0], we fall 

back to an inconsistency for the equations of motion. It can be seen in [25] that with this choice, an 
undesired relationship appears between the field equations, for the components 1 — 2(r — 0), which 
forces the theory to be restricted to TT, where /(T) is linear in T. So we make here another choice 
of tetrads. One possibility, consistent with the equations of motion, is that for which we obtain a set 
of tetrads through a Lorentz transformation for the diagonal case, leading to [26]. 


{e\} 


0 0 

0 sin 6 cos (j) r cos 6 cos (j) 

0 sin 6 sin cj) r cos 6 sin (j) 

0 e^^^cos0 —rsin0 


0 

—r sin0sin(/> 
r sin 6 cos (j) 
0 


(2.14) 


where we define the determinant of the tetrad by e = c?et[e“^] = sin0. With these tetrads 

we can reconstruct the metric with spherical symmetry 

dS'2 = [d02+sin2 ^ ('2.15) 


where the metric parameters {a(r),6(r)} are assumed to be functions of radial coordinate r and are 
not time dependent. 

Let us calculate the non-zero components of the torsion (2.4) and contorsion (2.5) 


rpU _ 

10 — 


9 o — 1 in ^ 

rfi2 _ _ Z;^1U _ 

^ 21 — ^ 31 ~ ^ ~ 


2 ’ zi oi r ’ ‘ 2 

We also calculate the non-zero components of the tensor in (2.6), giving 


99 nn e ~ 1) 

^22^ = ^ -V-; ^ (12 26) 


c 01 _ 

Dq — 


5-&(efe/2 1) 


J ^2 


C 12 _ Q 13 _ 
Oo — Oq - 


e"^ (a'r - + 2) 


4r 


(2.17) 


Through the torsion tensor components (2.16) and of (2.17), we can calculate the torsion scalar 
(2.7), which results in 



r 



(2.18) 


In general, the torsion scalar depends on the radial coordinate, as we will see below. 

Let us define the 4-potential electric being a 1-form A = A^dx^, where the components are 
given by {Afj^} = {^o(’'), 0, 0,0} and Aq is the electric potential because of the existence of just one 
charge. Let us define the Maxwell tensor as a 2-form F = A dx'^, where the components are 

anti-symmetric and given by F^i, = S/fiA^ — VvA^. With the specific choice of the spherical symmetry 
and the metric being static, there exists just an unique component of the Maxwell tensor Tjo (r) [27] . 
The field equations for using the Lagrangian density (2.8), are given by 


BCned 

= da 

Ij^a^dCNED 

[ dF \ 


[ dF \ 


whose solution, within ^ = 0, is given by 


(2.19) 


( 2 . 20 ) 


where q is the electric charge. 
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Through the above result, the field equations (2.12), with (2.13), are given by 


2^ (e''/2 _ t'/tt + ^ [b'r + - l) (aV + 2) 


fr + 7: = -‘^^NED - 2 — 
2 


(BCmedY. 


ej'’ - 
r 

e-^ r 


(.•/^ - 2) «V + 2 («‘/= -!)]/, + { = -2£„,„ - 2^ i?^) 


-1 


aV + 2 ( 1 - 


(l - e^/^) 


.-b r 


+53 -2”"+ 


+ {2b' + 4a'e''''^ - 6n') r - 4e'' + Se*'''^ - 4 


/r + -^ = —‘^•jC.ned ■ 


9F 

( 2 . 22 ) 


(2.23) 


In the next section we will re-obtain some of the most famous regular black holes solutions in 


GR. 


3 Recovering some known solutions in General Relativity 


In this section we will show that one can re-obtain all the regular black holes solutions in GR. To this 
end, we must demonstrate that the equations of motion (2.21)-(2.23) fall again to the same differential 
equations of GR. Then, we will stat by the particular case of the TT, where /(r) = T(r), frir) = 1 and 
fTrir) = 0. Hence, by subtracting the equation (2.22) from (2.21), one gets the following restriction 
a(r) = —h{r) + oq, with oq £ 5ft. From the symmetry of the temporal translation, one can redefine 
the temporal coordinate for obtaining this constant oq. Let us then set oq = 0, without loss of 
generality, such that b(r) = —a{r). As usually done for the case of regular black holes, let us rewrite 
the metric function a('r), in terms of another function depending on the radial coordinate as follows 


a(r) = —b{r) = In 


1 - 


2M{r) 


Now, by substituting all this into the field equations (2.21)-(2.23), one gets 


^NED 


( O Cned 
dF 


M'{r) 


= 0 . 


M"(r) „ 

‘2Cned H -= 0 . 


(3.1) 


(3.2) 

(3.3) 


These equations are exactly the same as that of GR, coupled with the NED for regular black holes 
(see the equations (7) and (8) of [28]^). Thus, all the solutions obtained in GR can be re-obtained 
through these equations here, as for example taking : M{r) = /{ 2 mr)^ which we will use later, 

nevertheless, in the context of the non-linear f{T) theory. 


4 Analogy between the known solutions and new black hole solutions 

The differential equations of motion for the f{T) theory coupled with the NED source, are differential 
non-linear coupled equations, rendering difficult how to obtain exact solution. But here, we will 
present a powerful algorithm for obtaining exact solutions for this theory. The basic idea is to use the 
Lagrangian density of NED, £ned, and it derivative with respect to the scalar E, dCNEo/dF = Cp, 
as functions that satisfying the field equations. We can impose that Cned satisfy (2.21), which by 
symmetry, and others considerations to be seen later, also satisfies (2.22). And finally, Cp has to 
satisfy the ultimate equation of motion (2.23). When the unique field component is given by (2.20), 
we will see that all the equations of motion are satisfied. But the expressions of Cppp and Cp, 
satisfying the set of equations (2.21)-(2.23), should not be related, one as derivative with respect to 

^Note that we use a definition of energy-momentum tensor in (2.13) different from the one used by Balart and 
Vagenas in the equation (6) of [28]. 
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F of the other. In such a situation, we impose this relation in order to obtain a consistent solution 
with this theory, leading to an equation with constraint, which has to be necessarily satisfied. 

Now we will present the algorithm. By subtracting the equation (2.22) from (2.21) one gets the 
following differential equation 


iv«MT)] + 


a'+ b' 
2(ef>/2 - 1) 


= 0 . 


By integrating this equation one gets 


frir) = exp 


1 r a'ir)+b'ir) \ 

2 J [eKr)/2 _ 1] ■ 


(4.1) 


(4.2) 


Note here that when we choose the relation b{r) = —a{r), from (4.2), one gets frir) = 1, which, 
by integration on T, leads to /(T) = T + To, where Tq G 5ft, known as the TT theory with cosmological 
constant, conform to the no-go theorem of [23]. Then, in order to do not recover the TT, we always 
have to ensure that b{r) ^ —a(r). Thus, we have here two possibilities to use the algorithm: a) fixing 
a functional form to the metric a(r), which is usual for regular solutions; b) fixing the functional form 
for the metric b(r). We will do it the coining section. 


4.1 Fixing a(r) 

In order to obtain new regular black holes solutions, it is common to rewrite the metric function a(r) 
as 


a(r) = In 



2M(r) 

r 


(4.3) 


where M{r) is an analytic function of the radial coordinate r, which falls to the black hole mass in 
the limit of spacial infinity, regarding asymptotically flat solution. Now we set (4.2) and (4.3) in the 
equation (2.21), for obtaining Cned, 


^NED = 


/W_2 ft 




2e-'’M/r(r) 


e^Fy^Mir) + r(l - M'(r)) j J / (r - 2M(r))) 


(4.4) 


where Cf = dCNEo/dF and M'{r) = dM{r)/dr. By subtracting (2.23) from (2.21), one gets Cf 


r2(2r2-2e'’«/V-2e'’«r2 
^2e^^yy'^r^ + ?,e^y^rM{r) - Se^^yy^rMir) - AM^{r) + Ge^y^M'^ir) - 

+8e^y^y^M^{r) + 2r^b'{r) - 2ey^y'^r^b'{r) - 6r^M{r)b'{r) + 7eyyyr'^M{r)b'{r) 

+ArM'^{r)b'{r) — 6ey^y^M^{r)b'{r) — 8r^M'{r) + 8e^*'’’^^^r^M'(r) -|- 8rM(r)M'{r) — 12e^yy^rM{r)M'(r) 
—2r^b'{r)M'(r) + e^yy^r^b'{r)M'{r) + M{r)b'{r)M'(r) — 2e^*'’'^/^r^M(r)6'(r)M(r) -|- 4r^(M'(r))^ 


Cf = — 


^q^e-yyfrir) - l) (r - 2M{r)f 


-2e’‘yy^r^{M'{r)f + 2r^M''{r) - 2ey^y\^M''{r) - 4r^M{r)M''{r) + 4e^(’')/V2M(r)M"(r)) 


(4.5) 


Now,, with the expressions (4.2), (4.3), (4.4) and (4.5) the set of the differential equations (2.21)- 
(2.23) are satisfied. We can obtain the non-null components of the Maxwell tensor, and taking (4.5) 
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in (2.20) one gets 


F^°{r) = |M2(r) 4 - g^b{r) _ g^3b{r)/2 ^ _ 2)rb'{r) 

-r^ [r6'(r) (2 - _ 2)M'(r)) + 2(4(e'’(’’)/2 - l)M'(r) - _ 2)M'\r) 

+ (eK’')/2 _ i)(_i + ghW _ ^M"(r)))] + rM(r) [r6'(r) (o - + 2(e^W/2 - 2)M'{r)'^ 

+4 (( 3 e '’(’')/2 _ 2)M'(r) + _ i)(2e''(’') - rM'{r))) ] I / [4g(e''(’')/2 _ i)(r - 2M(r))^('lj6) 


Now, we have to impose the relation between the Lagrangian density and its derivative in order 
to get a consistent solution. We then have the following consistent relation 


Cpir) 


dCNED f dF\ 
dr \dr) ’ 


(4.7) 


where F = (l/4)F''‘'F^i^. The imposition (4.7), within (4.4)-(4.6), yields the following differential 
equation 


[M(r) (2 — 2rb'{r)) + r {rb'{r) — 2M'{r))] ^1 — 26^*-'’^'^^^ M{r) + r + M'(r) ~ = 0(4.8) 


The equation (4.8) possesses two solutions 


b{r) = —a(r) , (4.9) 

and still, 

b(r) = 21n|e-“M 

with a(r) given in (4.3). 

Note now that there exist only two solutions which make the solutions of the differential equations 
being consistent, and satisfying (4.7). The first of these solutions, given by (4.9), falls into, where 
friT) = 1 in (4.2), because b{r) = —a{r), and then does not give any interesting thing, being just 
an analogous solution to that of GR. Let us now focus out attention to the second solution given in 
(4.10). 

We can now specify some models for M(r’), which is commonly done for regular black holes in 
GR [28]. Taking a simple model M{r) = me~'^ /{ 2 mr)^ with m € JR+, we can see that the black 
holes horizon is given by g^^irn) = = 0, with rp = —[q^/{2rn)]W[—q^/{2ra)'^], where W[y\ 

is the Lambert function [29], which solves lT[?/]e’^[*'l = x. We have a non-degenerate horizon for 
\q\y/e/2 < TO, then g^^{r) = 0 has two real and different roots; for to = \q\^l2 we have a degenerate 
horizon, and g^^(r) = 0 has a single real root; finally we have no horizon in the case |g|y^/2 > to and 
g^^ir) = 0 has no real root, we represent g^^(r) in figure 1. 

We still see that this is an asymptotically minkowskian solution, because in the limit of infinite 
space one gets = 1. The torsion scalar T(r) can be obtained from (2.18), yielding 

Qq^ F 8q^r (bm — 4e'^ + 8r^ (sm^ — 8e'* !jjyf _|_ 4g9 /mr^ 2 \ 

T{r) = - ^^^ -5- (4.11) 

-I- 2r (to — 

The torsion scalar (4.11) goes to zero in the limits r —0 and r -too. This shows that the torsion 
scalar is regular at the origem of the radial coordinate, also indicating a regular solution. 

As in some cases of charged solution, the torsion scalar can be regular, but it can exist non- 
detected singularity by this method [9], and we will calculate the curvature scalar and other ones 


M'{r) - 


M{r) 


- 1 


(4.10) 
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FigurG 1. Graphical representation of the function exp[—6(r)] for the cases \q\y/el"l < m (green), |g|y'^/2 = m (blue) and 
\q\^/e/2 > m (red). 


which exhibit singularity in GR. We can calculate the curvature scalar R = g^°‘ Rfiiyap, where 
R^ivap are the components of Riemann tensor calculated from the Levi-Civita connection^ [30, 36], 
by the expression R = —T — [30]. Thus, with the components (2.16) and the torsion scalar 

(4.11), we get the curvature scalar given by 


R{r) = 


4 ^2m — e'^ 4m?q^r + /(2™’')qr4 _|_ — rnq^{q^ + 4 e'* /( 2 ”n’),. 2 ^ 


mr [g 2 -g 2r{m — ' 


(4.12) 


The curvature scalar R{r) in (4.12) converges to zero for two limits r —>■ 0 and r —>■ +c», indicating 
that this solution is regular in all the spacetime. This can be confirmed by the Kretschmann scalar 
R^p^^ Qalculated also for the Levi-Civita connection. This scalar converge to zero for the two 
limits r —>■ 0 and r —>■ -l-oo. 

Here there is a huge algebraic subtlety to this class of solutions. First we note that the denom¬ 
inator of (4.11) and (4.12), as well as for the scalar R^iuapi vanishes with the same roots. So 

although the scalar being regular in the regions commonly analyzed, infinite space and origin, they 
diverge in interiot regions to the event horizon. We may check this by placing the curvature scalar, 
one can be made the same for the other scalars, only in terms of the function M{r) 


R{r) 


2[r — 2M{r)] — 2rM{r)[M'{r) + rM"(r)] -|- r^[M'{r)^ + rM"{r)]^ 

{M{r) + r[M'(r) — 1]}^ 


(4.13) 


with M{r) given specifically for each model, so it is true for this class of solutions in general. We can 
see in (4.13), this also occurs on (4.12), the denominator vanishes to a certain finite value 0 < Tg < rn- 
This shows us that there is a singularity in the region r = r^, which is located inside the event horizon, 
but before the origin. So even if all scalar are regular at the origin, this singular region (hipper-surface), 
prevents us from accessing it. The causal structure is then being identical to the Schwarzschild black 
hole in GR, showing that in fact this class solutions are black holes with a singular region covered by 
the event horizon. This same behavior occurs regardless of the model chosen for M{r). 

For this first class of solutions, the metrics functions (4.3) and (4.10) with M(r) = me~‘^ Rimr) ^ 
we represent functions exp[a(r)],exp[— 6 (r)], the curvature scalar R[r) and the Kretschmann scalar 
in Figure 2, for the non-degenerated case. We do the same in the figure 3 for the degen¬ 
erated case (extrem case). We note that although we have the restriction h(r) ^ —a{r), the zeros for 
exp[a(r)] and exp[— 6 ( 7 -)] still coincides, but the curvature scalar diverges in a value less than the event 

^The Levi-Civita connection is symmetric on (pi') indexes. The Riemann tensor calculated from 

the Weitzenbock connection in (2.3) vanishes identically [30, 36]. 
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FigurG 2. Graphical representation of functions exp[a(r‘)], exp[—6(r)] (left), and the functions R{r) ^ Kre{r) (right), in the 
non-degenerated case, for the values {m — q — 10}. 




FigurG 3. Graphical representation of functions exp[a(r‘)], exp[—6(r)] (left), and the functions R{r) ^ Kre{r) (right), in the 
degenerated case, for the values {m = y/eqj2^ q — 10}. 


horizon, characterizing a black hole. In Figure 2 we have the zeros of {exp[a(r)], exp[—6(r)]} on {r_ = 
2.322025582334978, r// = 13.98981145146388}, and the divergences (zeros for Kre~^{r)}) on 

{r^i = 1.4719210245449492, rs 2 = 8.906399770461864}. Then we found that the actual causal struc¬ 
ture is formed by an outer region of the event horizon r > rn, horizon himself r = rn, and other 
interior region rs 2 < r < rn with a singularity on r = rs 2 - All region r < rs 2 is excluded due to the 
impossibility extension of space-time. These are new class of black holes in solutions to /(T) Theory. 

The extreme case also results in new singular solutions. In figure 3 we see that there is a 
single zero to the functions {exp[a(r)],exp[—6(r)]} on VeH = 6.06530645576941, and a divergence in 
rs 3 = 2.505407627397168. There are a very high value, the order of 10^° to the Kretschmann scalar 
on horizon, but in the singularity, on r = r^a, the order is 10“^^, what appears to us that this is a 
solution with regular horizon. 

We calculate numerically the function frir) on (4.2), using (4.3) and (4.10) for this model, then 
we use the chain rule df/dT = {df / dr){dr / dT) (integrating to determine f{r) = f fxidT / dr)dr) 
obtain at the end the function /(r), which we can represent using a parametric graph in Figure 4. 
We can see in Figure 4 that the function f{T) does not depend linearly on T, then we are not address 
the TT here. 

The conclusion of this section is that the same behavior for this class of solutions can be shown 
in general, for any model for M (r), as listed in Table 1 in [28] . This shows that make a direct analogy 
between regular solutions to GR /(T) Gravity, results in singular black holes solutions here. This is 
due to strong restriction hir) ^ —air) arising from the No-Go Theorem [23], modifying the expression 
of the curvature scalar of the GR, as we see in 4.13. 


9 - 

























FigurG 4. Parametric representation to /(T) X T, in the non-degenerated case (left, r G [14, 500]) for the values {m — q — 10}, 
and degenerated case (right, r G [6, 50]), for the values {m = y/eqj2, q — 10}. 


4.2 Fixing b{r) 

We can now perform the same process of the previous section, but now, instead of fixing the functional 
form of a(r), we directly fix the functional form of b{r). This should be done with b{r) = — ln[l — 
(2M(r)/r)], but this leads to a complicated differential equation. Therefore, we will take into account 
the form 


b{r) = In 



2 M(r) 

r 


(4.14) 


Now we can in the same way as before, and in order to satisfy (2.21), with (4.2), we obtain 
^NEoir)- Afterword we get Cpir) in order to satisfy (2.23). Thus, with (2.20), we get the component 
from which, using the imposition (4.7), leads to the following differential equation 


AM (r) + r- I \jl -(2 + ra') — 2 


[2M(r)(l - ra') + r{ra' - 2M'(r))] = 0 . (4.15) 


Then, we have two solutions for (4.15), where the first is a{r) = —b{r), which, from (4.2), leads to 
TT, with friT) = 1. This first solution is not interesting here, because of its analogy with the GR. 
The second solution of (4.15) is given by 


= -2 / J L 


eKG/2 _ x' 


dr . 


(4.16) 


As we have to obtain regular black hole solution, the horizon is defined by g^^irn) = e = q, 

which is well described for 


M(r) = — 


r — 2m(r) 


(4.17) 


because, for this, one has e = 1 — [2m(r)/r]. 

First we analyse in general the scalar of Kretschmann for this class solutions 



-30 + 41 


r — 2m(r) 


+3m'{r) + 2 - 2^;-^) ] } ■ 
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(4.18) 





























We note that solutions can only be analytic throughout the space if the restriction r — 2m(r) > 0 is 
satisfied. As cases of non-degenerate horizons break this condition, we will focus our attention only 
to the extreme case where this condition is valid. 

Let us take again one model for mlr). The model for m{r) is given by m(r) = moe“'^ /( 2 moi’)^ 
with mo € 5J+. The events horizon of this solution is given by th = —[q^/ (2mo)]W[—(2mo)^], with 
This solution is asymptotically flat and regular in the spacial infinity. The curvature 
scalar is given by 


R(r) = - 


l(2mor) 


2 r 


mo — e'^ ^ g(j2/(2mor)j, ^ — 2mo + 


.(4.19) 


The curvature scalar and R^vap tend to zero in the limit r —>■ 0, showing that the solution is 

regular at the origin of the radial coordinate. 

Let’s look again at the metrics and scalar functions, which represent the figure 5. We can see 
that we have only one horizon rf^H 2 = 30.326533058652203. The torsion scalar r(r), the curvature 
scalar R{r) and the Kretschmann scalar They are finite in all space-time. 

We again calculate numerically the function /t(?') on (4.2), using (4.16) and (4.14) for this model, 
then we use the chain rule df /dT = {df /dr){dr/dT) (integrating to determine /(r) = J fxidT/dr)dr) 
obtain at the end the function /(r), which we can represent using a parametric graph in Figure 6. 
We can see in Figure 6 that the function f{T) does not depend linearly on T, then we are not address 
the TT here. 


5 Conclusion 

As well as in GR, by analogy, we establish the equations of motion for the /(T) theory with matter 
source of non-linear electrodynamics (NED), with spherical symmetry in 4D, in order to obtain regular 
black holes solutions. Through the equations of motion we show that for the particular case where 
/(T) = T, we can re-obtain all the solutions previously obtained in GR. Also, through the previous 
the equations of motion, we show that the unique regular black holes solutions which can be obtained 
in /(T) theory, with the non-linear terms in the function and the tetrads chosen, are in conformity 
with the no-go theorem in [23] for h{r) ^ —a{r). 

We elaborate here a powerful method for obtaining new regular black holes solutions, in the 
case of /(T) theory with the NED source. Thus, by fixing first the functional form of the metric 
a(r), we get a new class of black hole (singular), where the scalars, torsion scalar, curvature scalar 
and Rfivap, converge to zero at the limit of the origin of the radial coordinate, but diverge 

for the value rs2 < rn- In this case we choose just one model where M{r) = me~'^ /{2mr)^ with 
a{r) = ln[l — (2M(r)/r)] and b{r) = 21n{[M'(r) — {M{r)/r) — a(r)]/a(r)}. In this case the solution 
is asymptotically flat, possesses an event horizon covering a singular region r = rs 2 - Still study the 
extreme case, where m = y/eql 2 , resulting in the same structure. Both cases present the same causal 
structure of the Schwarzschild black hole on GR. 

We also fix the functional form of b{r) = — ln[l —(2m(r’)/r)], obtaining a new class of regular black 
hole solution. The solution with m{r) = moe~'^ /(2mo’') jg asymptotically flat and regular everywhere. 
They possess a degenerate horizon behaving as the ones corresponding to GR, but here b{r) ^ —ci{r) 
and a different curvature scalar R{r) given by (4.19). 

Our perspectives about these new class of charged regular solutions that it is still possible to 
find another solutions here. We also believe these classes of charged regular solutions must obey the 
energy conditions for f{T) Gravity [32], but breaks specifically the SEG, as well as in GR [33] (see 
also [34, 35] where the DEG and WEC are satisfied). We also think there should be some means to 
test the requirement that electrically charged solutions do not have the correct weak field limit as in 
GR [37]. This should be thoroughly studied on another paper. 
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